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Abstract. We investigate the properties of KOV stars 
with Hipparcos parallaxes and spectral types taken from 
the Michigan Spectral Survey. The sample of 200 objects 
allows the empirical investigation of the magnitude selec- 
tion (Malmquist) bias, which appears clearly present. By 
selecting those objects that are not affected by bias, we 
find a mean absolute magnitude of M v — 5.7, a downward 
revision from 5.9 mag. listed in Schmidt-Kaler (1982). 
Some objects have absolute magnitudes far brighter than 
M v — 5.7, and it is suggested that these objects (w 20% 
of the total sample) are KOIV stars which may have 
been mis-classified as a KOV star. The presence of the 
Malmquist bias in even this high quality sample suggests 
that no sample can be expected to be bias-free. 



Key words: Stars: distances - 
eters - Stars: late-type 



Stars: fundamental param- 



1. Introduction 

The Hipparcos trigonometric parallax measurements of 
more than 100 000 stars (ESA, 1997) provide an excellent 
basis to determine the fundamental parameters of stars. 
Yet, some, not always trivial, problems arise which make 
the conversion from the measured parallax to intrinsic ab- 
solute magnitude of an object not straightforward. For 
example, the Lutz-Kelker effect (Lutz & Kelker, 1973), 
results in too faint magnitudes for large relative errors 
<j/tt, while the Malmquist bias results in too bright mean 
absolute magnitudes, because at the observed magnitude 
limit, brighter objects will be included in a sample, while 
fainter objects will not. Additional complications are listed 
in Brown et al. (1997). 

In a previous paper (Oudmaijer, Groenewegen & 
Schrijver 1998 - hereafter OGS98) we have shown em- 
pirically that the Lutz-Kelker bias is present in trigono- 
metric parallaxes. This was done by comparing the best 
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Hipparcos parallaxes {a /it < 5% - defining a 'true' par- 
allax sample) with lower quality ground-based parallaxes 
of a large sample of stars. The data showed that, for in- 
creasing ct/tt, the derived absolute magnitude of an object 
indeed becomes too faint in a manner consistent with the 
Lutz-Kelker predictions (see also Koen 1992), but for even 
larger a/n, the derived magnitudes became too bright by 
up to 2 magnitudes. The sample was evidently not ham- 
pered by only one type of bias, but by at least two. The 
first being the Lutz-Kelker bias, the second was called the 
'completeness effect', which we now identify as the mag- 
nitude selection Malmquist bias. 

To investigate this further, we tackle the problem in 
a similar way by analyzing a sample of stars for which it 
may be hoped that all have the same intrinsic magnitude. 
To this end, we have investigated a sample of stars with 
well-defined spectral types, the KOV stars. 

2. Sample selection 

To determine the absolute magnitudes of stars with the 
same spectral type, a coherent and homogeneous database 
of spectral types is needed. The Michigan Spectral Survey 
Volumes 1..4 (MSS, Houk & Cowley 1975; Houk 1978; 
Houk 1982 and Houk & Smith-Moore 1988, providing 
spectral types of the HD catalogue in Declination from 
-90° to -12°) is such a database. We chose to investigate 
KOV stars, as these objects are relatively close by and 
will not suffer much from interstellar extinction, while the 
number of objects is relatively large. The selection criteria 
from the Hipparcos Catalogue (ESA, 1997) were: 

(i) Spectral type = 'KOV ' (Field H76, the sources for the 
spectral types are listed in Field H77) 

(ii) Goodness-of-fit < 3 (Field H29) 

(hi) Number of rejected data < 10% (Field H30) 

563 objects in the Hipparcos catalogue have 'KOV ' 
listed in their spectral type entry, but the majority of the 
stars has the spectral type taken from other sources than 
the MSS or are listed as 'G8V/K0V. These objects were 
rejected, leaving a sample of 201 objects. Only one of these 
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Fig. 1. Derived M v as function of several parameters. Er- 
rorbars on the absolute magnitudes are for convenience 
only shown in the lower panel, and are often smaller than 
the plotsymbols. The solid lines are drawn according to 
Eq. 1 in OGS98, and explained in the text. 

has a negative parallax (HD 219882), and was also rejected 
for further analysis, one object (HD 170132) has no (B-V) 
listed in the Hipparcos catalogue, its value was taken from 
the simbad database. 

The selection thus yielded 200 objects. The average 
error on the parallax and its scatter, are 1.4 ± 0.6 mas, and 
the bulk of the sample has parallaxes larger than 10 mas, 
probing the nearest 100 pc. The quality of the parallaxes 
is extremely high; the average a/n is 11%, indicating a 9cr 
detection on average. 

3. Properties of the sample 

The absolute magnitude, derived from the parallax and 
the F magnitude, neglecting interstellar extinction, is plot- 
ted in Fig. 0. The upper panel shows M v against (B-V). 
The unweighted mean M v = 5.06 ± 1.26 (the r.m.s. de- 
viation around the mean) is almost 1 magnitude brighter 
compared to what is expected for KOV stars (M v =5.9, 
Schmidt-Kaler 1982, hereafter SK82). Some objects are 
even 6-8 magnitudes brighter than a normal KOV star. 
A trend in (B-V) may be present, as the redder stars cor- 
respond to the intrinsically brightest objects. The relation 
between M v and V (middle panel) shows a large scat- 
ter which seems to increase for fainter objects. There is a 
strong correlation between the derived M v and the mea- 
sured parallax (lower panel). For small parallaxes, the in- 
trinsic magnitude becomes brighter and, interestingly, for 
the smallest parallaxes, no objects have intrinsic magni- 
tudes that are even close to M v = 5.9. A strong limit to 
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Fig. 2. The distribution of M v , binned to 0.25 mag. bins. 
The dotted line represents a two component Gaussian fit 
to the data. 

the derived M v as function of parallax is present, which 
is due to the 'completeness effect' mentioned in OGS98. 

As discussed in OGS98, the difference between the de- 
rived absolute magnitude of an object from its parallax 
and the limiting observed V magnitude of a sample define 
a 'forbidden' region. The reason is that stars that would 
have been present in the fainter regions (the lower left 
hand corner of the lower panel) are simply too faint to be 
included in the sample. An upper bound is also present, 
reflecting the fact that fewer bright objects are present 
than faint objects. The solid lines in the figure indicate 
the regions where no data are expected, and are drawn 
according to Eq. 1 in OGS98, with limiting magnitudes 
corresponding to the faintest V magnitude in this sample, 
V = 10.19, respectively the brightest, V = 5.95. 

We now identify this 'completeness effect' with the 
magnitude selection Malmquist bias, which is working on 
exactly the same principle, and effectively forbids the use 
of entire samples to derive their mean absolute magni- 
tude, without correcting for it and/or investigating when 
the bias starts to dominate. It has been advocated to first 
plot luminosities as a function of distance, a parameter 
directly related to the distance such as the red-shift in 
the case of galaxies (e.g. Sandage 1994) or, in this case 
parallax, to assess the presence of selection biases. Such 
diagnostic plots, sometimes called Spaenhauer diagrams 
after Spaenhauer (1978), also serve to identify Lutz-Kelker 
type biases (see OGS98). Sandage (1994) showed that his 
sample of galaxies suffers from the Malmquist bias (with 
AM = 1.386 x a 2 , for a uniform space distribution, with 
a the assumed intrinsic scatter of the absolute magnitude 
distribution, see e.g. Hanson, 1979) when he compared 
the average with a sub-sample, easily identifiable in the 
diagrams, which is not affected. 

Let us now derive the Malmquist correction for our 
sample, Fig. [l] shows that those objects with n > 20 mas 
are not affected - the unweighted mean of these objects 
returns a value of M v = 5.69 with an (intrinsic) scat- 
ter of 0.4 mag. The entire sample returns an average of 
5.06 ± 1.26 mag. The expected Malmquist correction is 
0.22 mag for a — 0.4, so the difference between the de- 
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Fig. 3. The resulting KOV sample - 159 objects with 
M v > 4.5. The solid lines are as in the previous figure. The 
triangles indicate the mean and its scatter in the parallax 
bins 10-12.5; 12.5-15; 15-17.5; 17.5-20; 20-25; 25-30; 30-35; 
35-60 and > 60 mas. 
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Fig. 4. Cumulative parallax distribution of the KOV sam- 
ple. The solid line indicates a fit through the data points 
with tt > 20 mas. The bins are 0.05 wide in log(7r) units, 
the errorbars represent the statistical error (VW). 



rived M v for the entire sample and that of the unaffected 
sub-sample is much larger than what the Malmquist bias 
predicts. This is at first sight hard to understand, but may 
be related to the question why we would find KOV stars 
which are up to 6-8 magnitudes brighter than expected. 
Apart from the rather unlikely possibilities that parallax 
errors would result in such deviant values (these objects 
have very high signal-to-noise detections) or that the class 
of KOV stars can have such a large range of intrinsic mag- 
nitudes, it may be more likely that the discrepancy is due 
to spectral misclassification. 

Some information may be gained from Fig. |^, where 
the distribution of the derived M v values is shown. The 
distribution peaks close to 5.7, but is not symmetric 
around the mean; a secondary maximum appears close 
to M v =3.5. The presence of a secondary peak strongly 
suggests that an additional population of stars is present. 
These could be objects with a different spectral type as 
the peak roughly agrees with the magnitudes of K0IV stars 
(M v = 3.1, SK82), while K0III giants may also be present 
(M v = 0.7, SK82). It is hard to make a good distinction 
between the different groups, judging the gaps between 
the M v w 5.7 and M v w 3.5 objects in Figs. [I] and ||, the 
separation seems to be present for M v = 4.5. 

159 of the 200 objects are present in the 'faint' sam- 
ple with My > 4.5, the remaining 41 stars have brighter 
intrinsic magnitudes. The average (B-V) of the faint (i.e. 
KOV) sample is 0.82 ± 0.05, consistent with the intrinsic 
colours for the group (0.81, SK82), also suggesting that 
our neglect of interstellar reddening is warranted. The re- 
maining objects have a redder average (B-V) of 0.94 but 
with a large scatter of 0.16 mag. If we reject the 6 bright- 
est objects in this sample, the scatter is reduced and the 
average (B-V) = 0.89 ± 0.08, with an average M v of 
3.5 ± 0.5, consistent with a K0IV nature of the sample. 
Although Schmidt-Kaler (1982) does not list the (B-V)o 
for K0IV objects, the interpolated value between KOV and 
K0III stars is 0.90, close to what is measured. The remain- 
ing 6 objects have even redder colours, (B-V) = 1.2 ± 0.2 



with an average M v — 0.1 ± 1 mag, suggesting that these 
may be K0III stars. 

The simplest explanation for the large range in ab- 
solute magnitudes then appears that the sample of KOV 
stars in the MSS survey is contaminated by K0IV objects 
(about 20% of the entire sample), and perhaps suffers from 
contamination from intrinsically even brighter objects. 

4. The KOV subsample and the Malmquist bias 

In the following, we will continue with the KOV sample, i.e. 
the 159 objects with inferred intrinsic magnitudes fainter 
than 4.5. Our interest is whether the 'completeness effect' 
or Malmquist bias affects the determination of the intrin- 
sic magnitude of the sample under consideration. Fig. |^ 
shows the same figure as the lower panel of Fig. |], but now 
only for the KOV sample. As before, there is a clear trend 
visible. The smaller the parallax, the brighter the mean 
is. We have binned the data in steps of tt, and calculated 
the mean and its scatter. In the interval tt = 10-12.5 mas 
(which are still 6-lOc detections), the mean is 0.7 mag 
brighter than in the interval 60-100 mas. 

The change in mean absolute magnitude is easily un- 
derstood. This can be learned from the volume complete- 
ness of the sample. Fig ^| shows the cumulative parallax 
distribution of the stars. The distribution flattens below 
tt < 20 mas indicating that the sample is complete to 20 
mas. The solid line represents a weighted least-squares fit 
to the data between 20 and 100 mas, with slope — 3±0.15, 
implying a uniform space distribution, consistent with 
the small volume probed to 50 pc. The Malmquist bias 
only occurs for volume-incomplete samples, and indeed, 
for tt > 20 mas, the mean magnitude in the bins does not 
change in Fig. ||, but it does for the lower parallax values. 

What is the effect in this particular case on the de- 
rived absolute magnitudes of KOV stars if the Malmquist 
bias would not have been taken into account? The un- 
weighted mean of the 159 objects is 5.59 ± 0.42 mag, while 
the unweighted mean for the unaffected sample (85 stars 
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with 7T > 20 mas) is M v = 5.69 ± 0.40 mag. The scatter 
of w 0.40 reflects the intrinsic scatter rather than errors 
arising from the measurement uncertainties, and may for 
example be due to variations in metallicity, rotation pe- 
riod or unseen binaries. The difference between the two 
values is more in agreement with the prediction that the 
Malmquist bias is of order 0.2 mag - this is dependent on 
the distinction between the KOV and the 'K0IV samples, 
because a fainter cut-off value results in a slightly smaller 
scatter around the mean. 

However, there is one significant difference with e.g. 
the situation of red-shifts as distance determinations: in 
the parallax case the relative observational error er/7r, 
is much larger than in the case of the red-shifts, so a 
straightforward averaging of the derived absolute mag- 
nitudes should be replaced by a weighting scheme. This 
will decrease the effect of the Malmquist bias somewhat: 
The objects that are more prone to the selection effects 
are further away, and have larger relative errors on the 
parallax, they will therefore have less weight. Since the 
error is asymmetric in magnitudes, we now have to cal- 
culate the weighted mean in 'reduced parallax' (10°' 2Mv ). 
The weighting of all 159 objects now results in a mean 
M v = 5.65, while the 85 objects with 7r > 20 mas have a 
weighted mean of 5.69 mag. So, for this sample, consisting 
of both high quality parallax measurements and spectral 
types, not taking into account the Malmquist bias would 
result in an unweighted mean too bright consistent with 
the expected value of the Malmquist bias, and a weighted 
mean absolute magnitude that is too bright by 0.04 mag. 

The main result concerning the 'true' intrinsic mag- 
nitude of KOV stars is that the sample which is not af- 
fected by contamination by K0IV stars yields a value 0.2 
mag brighter than has been listed in the literature so far 
(SK82). The existing calibrations apparently need a re- 
vision and this work illustrates the power of Hipparcos 
trigonometric parallaxes. An additional result is that the 
stars that we tentatively identify as K0IV objects seem 
to be 0.4 mag. fainter than expected. However, we do not 
put much weight to this result, as these by implication 
would be mis-classified KOV stars, and thus likely to be 
those K0IV objects that are on the fainter side of the dis- 
tribution in the first place. A detailed study of this effect 
is beyond the scope of the present paper. 

5. Final remarks 

We have studied a sample of 200 KOV stars, taken from 
the best collection of spectral types available, the Michi- 
gan Spectral Survey, which have excellent trigonometric 
parallaxes from the Hipparcos mission. In our high qual- 
ity data, the Malmquist bias occurs already when the 
measured parallax is 20 mas, the point where volume- 
incompleteness sets in. The presence of the Malmquist 
bias is readily seen, when the Spaenhauer diagram is used 
as a diagnostic tool. Only because KOV stars are intrin- 



sically faint, the Malmquist bias starts to play a role this 
quickly. For intrinsically brighter objects, or samples with 
fainter limiting V magnitudes, the completeness limit will 
be pushed to lower observed parallaxes. 

So far, we have not discussed the Lutz-Kelker (1973) 
effect as shown in OGS98, this effect becomes impor- 
tant when the relative error on the parallax is still com- 
paritively low (10-20%). Since the absolute errors on the 
parallaxes are almost all the same (« 1.4 mas in our sam- 
ple), this corresponds to an observed ir of 7-14 mas. In the 
case of the faint KOV stars, the Malmquist effect domi- 
nates at these parallaxes, so objects that otherwise would 
have had too faint derived absolute magnitudes are ex- 
cluded from the sample, perhaps lowering the effect of the 
Malmquist bias. To assess such effects, one has to exam- 
ine intrinsically brighter objects, as done in OGS98, or 
by Kaltcheva & Knude (1998) who investigated B stars. 
The latter authors showed that for well determined par- 
allaxes (cr/n < 10%), the absolute magnitudes of B stars 
derived from the parallax agree with the absolute magni- 
tudes derived from R(3 photometric distances, but for less 
well determined parallaxes (<j/tt between 10% and 20%), 
the individual objects are too faint, in agreement with the 
prediction for Lutz-Kclker bias. 

As a final comment, it is expected that each sample 
of stars will prove to be sensitive to the Malmquist and 
Lutz-Kelker biases in its own unique way, and only careful 
examination of the data can make their, sometimes hid- 
den, effects visible. An additional result of this work is 
that around 41 out of 200 KOV stars may be misclassificd 
K0IV stars. 
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